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The experimental observation of lepton-number violating processes would unambiguously indicate the Majo- 
rana nature of neutrinos. Various AL = 2 processes for pseudoscalar meson M[ decays to another pseudoscalar 
meson M2 and two charged leptons Ci, ti {hT[ — > C^i^M^) have been studied extensively. Extending the exist- 
ing literature on the studies of these kinds of processes, we consider the rare decays of heavy mesons to a vector 
meson or a pseudoscalar meson. These processes have not been searched for experimentally, while they may 
have sizable decay rates. We calculate their branching fractions and propose to search for these decay modes in 
■ the current and forthcoming experiments, in particular at the LHCb. 

(N 

^ . 

I. INTRODUCTION 

rS 1' The neutrino oscillation experiments have proved that neutrinos are massive However, the nature of neutrino masses is 

' still one of the main puzzles in contemporary particle physics, i.e., are neutrinos Dirac or Majorana particles? As we all known 
that the Majorana mass term violates lepton number by two units (AL = 2). Thus, the unambiguous answer to the question above 
(— I is the experimental observation of a lepton-number violating (LV) process. 

I— ' , Various AL - 2 processes have been studied in the literature t5l-|Ml- Among them Atre et at. have studied 36 LV processes 
' from K, D, D,, and B decays, generically written by: 

O ■ < ^14^2' (1) 

t — ! 

(y-j , where Mf and £f (/ = 1 , 2) denote charged pseudoscalar mesons and leptons, respectively. 

' Most of these processes have been searched for and the non-observation in the current experiments set the bounds on branching 
fractions. In turn, they led to some stringent constrains on the mixing parameters between Majorana neutrino and charged lepton 
( <~>) directly. However, there are still more LV heavy meson decay modes that have not been studied experimentally, that may 
have sizable branching fractions in theory. In particular, heavy mesons (with c, b flavors) are easier to identify and the LHCb 
I , experiments will provide us with a large data sample. So as an extension of current existing calculations, we explore some new 
' AL = 2 decay modes in this paper We mainly consider the rare decays of heavy mesons D, D,, and B to vector meson final 
states. Since the LV heavy meson decay modes under our consideration have no experimental results, we cannot extract the 
mixing parameters through these decay channels like Atre et at. did. 
5^ \ However, those processes considered by Ali jsj] and Atre ||6|] are clearly correlated with the decay modes under our consider- 
ation, with the same mixing parameters specified by the charged lepton flavors. We thus adopt the numerical values of mixing 
parameters extracted from Ref. and the decay widths and branching fractions of heavy mesons for our processes can be 
predicted correspondingly. We choose the strongest constrains on mixing parameters from Ref. |j6t] as input in our study, in order 
to be conservative. 

we mainly consider the heavy pseudoscalar meson D, D^, and B to vector meson final states. From theoretical point of view, 
the decays of vector mesons may have different and uncoiTelated rates from that of the pseudoscalars if there is other type of 
new physics, like a heavy particle exchange of either a pseudoscalar/scalar or a vector boson. Thus, it is well motivated to carry 
out the complementary searches for all of the availabel final states. 

The paper is organized as follows. In Sec. II we outline the useful formulas to set the general notation. We list the constraints 
on mixing parameters and give the Monte Carlo sampling of the branching fractions as a function of the heavy neutrino mass in 
Sec. Ill, and draw our conclusion in Sec. IV 
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FIG. 1: Feynman diagram corresponding to the AL = 2 processes via Majorana neutrino exchange. 



n. THE GENERAL FORMALISM FOR LEPTON-NUMBER VIOLATING DECAY 



The simplest renormalizable extension of the standard model (SM) to generate neutrino Majorana masses is to introduce n 
right-handed SM singlet neutrinos NbR {b = 1,2, ■ ■ ■,«). Therefore, the complete neutrino mass sector is composed of both 
Dirac masses that produced via the Yukawa couplings to the Higgs doublet in the SM, and possible heavy Majorana mass term 

3 i: Nf^Bbb'Nb'R + h.c.. 

b,b'=l 

In terms of the mass eigenstates, the gauge interaction lagrangian of the charged currents now has the following form: 

r 3 T 3+/I 



£ = — —W^ 

V2 ^' 



V f=e m=l f=e m'=4 



+ h.c. 



(2) 



where Pi = ^(1 - 75), v„(m = 1, 2, 3) and N,„'{m' = 4, ■ ■ ■ , 3 + n) are the mass eigenstates, Ucm is the mixing matrix between 
the light flavor and light neutrinos, and Ve,,,' is the mixing matrix between the light flavor and heavy neutrinos. 

The basic process with AL = 2 shown in Fig. [1] with exchange of two virtual SM W bosons can be generically expressed by: 



^1 ^2 



(3) 



where {1^2 = e,fi, t. The process can occur only if neutrinos are Majorana particles. Unfortunately, the transition rate of this 
AL = 2 process encounters a severe suppression either due to the small neutrino mass like mli^JM^, or due to the small mixing 

I |2 

\V{,m' Vtim' I ■ However, when the heavy neutrino mass is kinematically accessible, the process may undergo a resonant production 
of the heavy neutrino, thus substantially enhancing the transition rate. The resonant contributions of heavy Majorana neutrinos 
to AL - 2 processes involving two charged leptons and another pseudoscalar meson have been considered in the Ref. |6]. In 
this paper, we extend the study of the heavy Majorana neutrinos to AL - 2 processes involving two charged leptons and a vector 
meson final state. 

The Feynman diagram for the LV decay of heavy meson Mi into two charged leptons {\, €2 and another meson M2: 



M\(qi) ^ (\{pi) q{p2) M-^{q2) . 



(4) 



is shown in Fig. |2] According to narrow-width approximation, the tree level decay amplitude when M2 is a pseudoscalar meson 
is given as [|6(]: 



Mf, til iiPRVe. 



Xqi - piY -ml + iTM^m4 
when M2 is a vector meson, the tree level decay amplitude can be written as J^: 



(qi - Pip- -m\ + iTN,m^ 



+ {p\ ^ Pi), 



+ (P\ ^ Pi), 



(5) 



(6) 



where Gf is Fermi constant; Vj^ is the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements; Z^, is the decay constant for 
meson M,-; q\,qi, p\, pi are the momenta of mesons M\ , M2 and leptons {\, I2, respectively. Here we consider the case when 
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FIG. 2: Feynman diagram corresponding to the lepton-number violating decays M\{qi) — > ^^CpO^J (^2)^2 (^2). 



TABLE I: The decay modes of the leptonic-number violating decays M^(q\) — > C\(p\)l^{p2)M2 (^2) and the ranges of the mixing parameters 
|Vf,4Vf,4|. The rages are calculated base on the correlative decay modes from the Ref. [60. 



Decay mode 


Mixing element |Vf,4Vfj4| 


Range of m4(MeV) 


Range of |Vf|4Vf,4| 


D+ -» e+e+p- 
^ e+e+if- 

B+ -> e+e+D*- 




776 - 1869 
892 - 1869 
776 - 1968 
892 - 1968 
1870-5278 
1969 - 5278 
2011 -5278 
2113 -5278 


0.002 - 0.998 
0.002 - 0.998 
0.023 - 0.992 
0.022 - 0.992 
0.236 - 0.992 
0.231 -0.992 
0.231 -0.992 
0.232 - 0.992 


B* iA*iA+D- 
B* ti*ti^D- 
B* fi+fi+D"- 

B* ^ yU+yU+D^ 




882 - 1863 
1975 -5173 
2074-5173 
2116-5173 
2218-5173 


0.002 - 0.993 
0.28 -0.971 
0.28 -0.971 
0.28 -0.971 

0.327 -0.971 


Z)+ e+yu+p" 

-> e*n*K'- 
D+ e+yu+p- 
Z)t ^ e+fi+K'- 

S+ e*iJ*D'- 
B+ e*n*D'- 


|VmV^4| 


776 - 1869 
892 - 1869 
776 - 1868 
892 - 1968 
1870 - 5278 
1969 - 5278 
2011 -5278 
2113 -5278 


0.02 - 0.497 
0.023 - 0.497 

0.02 - 0.497 
0.023 - 0.497 
0.242 - 0.49 
0.236 - 0.49 
0.236 - 0.49 
0.236 - 0.49 



only one heavy Majorana neutrino is kinematically accessible and denote it by A^^, with the corresponding mass m4 and mixing 
with charged lepton flavors V(n. is the total decay width of the heavy Majorana neutrino, summing over all accessible final 
states. Then the partial decay width F^' and the normalized branching fraction BR - F^' /Fm, for the LV process Eq. (|4|i can be 
calculated by the decay amplitude. Following the approach of Ref. la], we will take the mixing parameter Vca and the mass nin 
as phenomenological parameters. 

III. MONTE CARLO SAMPLING FOR LEPTON-NUMBER VIOLATING DECAYS 

The key step to calculate decay widths and branching fractions of the LV heavy meson decays is to determine the limits on the 
mixing parameters |yf,4yf24| and neutrino mass nin in Eq. (|5]l and Eq. (|6]l. Generally speaking, one can determine limits on the 
mixing parameters from the LV heavy meson decay modes which have the current experimental limits on branching fractions 
and determine the mass of neutrino by kinematics. However, as mentioned in the introduction, since the LV decay modes which 
we studied with vector meson and several pseudoscalar meson final states are missing in directly experimental searches, so we 
cannot yet get information of mixing parameters |yf,4yf24| from those decays. We thus propose the direct searches for those 
modes in the existing and forth coming experiments such as in CLEO, Z?-Factories, and the LHCb. On the other hand, there 
are direct experimental results on the processes that may share common mixing parameters with those under our consideration. 
These decay modes given by Ref. ^ and Ref. iH have been summarized and translated into the direct bounds in Ref. la] . 
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m^CGeV) 

FIG. 3: Theoretically excluded regions inside the curve for the branching fraction of e^e*M^ modes versus Majorana neutrino mass 

m^. Regions below the curve are theoretically allowed. The curves of £), — > eep, D — > eep, — > eeK' and D — > eeK', denoted by solid line, 
break line, dotted line and dash-dotted line, respectively. 

We thus adopt them for consistency and carry out our analyses. We first obtain the consistent limits on mixing parameters from 
some decay modes which have the current experimental bounds on the branching fractions for the heavy mesons with Ref. f^. 
We then translated the limits on mixing parameters to the relevant decays modes of the heavy mesons D, and B in the Table I. 
Depending on the flavors of the final state leptons, the mixing parameters probed are \Ve4[^, |V^4|^ and |y(.4V^4|, corresponding to 
the decay modes; 

M+ ^ e+e+Mj , ^ e+yU+Mj and ^ yU^+Mj , (7) 

respectively. We list the most stringent limits on IV^P, |V^4p and |ye4V'^4| for the 21 new decay modes in Table I. The ranges of 
heavy neutrino mass m4 in the Table I are determined by the kinematics accessible. 

When performing the calculations, the input pa rameters for the CKM matrix elements and the decay constants of pseudoscalar 
and vector mesons are chosen as follows iflOl [l2l - [l5ll : 

\Vub\ = 0.00359, \Vcd\ = 0.2256, |y,.,| = 0.97334, 

/d. = 0.2226 GeV, fof = 0.266 GeV, /b± = 0.190 GeV, 

/p± = 0.220 GeV, fx-* = 0.217 GeV, fo- = 0.31 GeV, fo-- = 0.315 GeV. (8) 

We note that there may be some eiTors in determining the decay constants ifl^ , but they would not result in any qualitative 
difference for our predictions for the SM-forbbiten modes. 

With these parameters and the limits on mixing parameters and coiTesponding mass ranges in the Table I, the decay widths 
and branching fractions of the heavy mesons D, and B are calculated coiTespondingly. We perform a Monte Carlo sampling 
of the branching fractions and the mass of the heavy neutrino, i.e., we plot the excluded region of the branching fractions as a 
function of m^, as shown in Fig.[3]~ Fig. |5] for the modes in Eq. O. The regions inside the curve are excluded by the direct 
experimental searches for the various LV decay modes of heavy mesons as obtained in Ref. |6]. The theoretical allowed regions 
are below the curve, i.e., the regions below the curve are the currently allowed branching fractions for those LV heavy meson 
decay modes in the Table I. 

From the figures, one can see that, if the heavy neutrino mass is located in the range from 1 GeV to 2 GeV, even with the 
most stringent constraints on mixing parameters, our theoretical predictions of upper bound of the branching fractions can be 
large, for example, Br{Ds -> eep) = 1.8 x lO"-*, Br{Ds — > e^ip) - 2.0 x lO"-' and Br{D -> e^p) = 1.3 x 10""*, since about 
2.4 X 10^ D^D 1 17] events and 5.5 x 10^ D*s^D^^ ifisll events have been collected by CLEO collaboration, these mentioned decay 
modes can be analyzed in the current experiment, which will provide us a strong information of 1 GeV to 2 GeV neutrino mass. 
But if the heavy neutrino mass is heavier, around 2 GeV to 4 GeV, we cannot use the D or decay modes to detect the LV 
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FIG. 4: Theoretically excluded regions inside the curve for the branching fraction of — > e^^* modes versus Majorana neutrino mass 
m^. Regions below the curve are theoretically allowed. The curves of — » ejip, D — > epp, £), — » e^K' and D — > e^K', denoted by solid line, 
break line, dotted line and dash-dotted line, respectively. 




FIG. 5: Theoretically excluded regions inside the curve for the branching fraction of — > ^^fi* modes versus Majorana neutrino mass 
m^. Regions below the curve are theoretically allowed. 



processes because of the kinetic limited, and the B decay modes are favored, whereas our predicted branching fractions of the 
LV B decay channels are lower than 8.0 x 10"^, which cannot be detected by current B-factories, but in the forthcoming LHC 
experiments, around 10" ~ 10'^ B meson events are expected [19] , all the LV B decays modes we studied can be effectively 
searched for by LHCb. 

As a final remark, we would like to reiterate the advantage of our treatment in searching for the resonant production and 
decay of a Majorana neutrino in vector meson decay. Although kinematically limited, the signal rate for the rare meson decay is 
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substantially enhanced due to the resonant nature. In contrast to the similar decay channels as discussed in Ref. f3\, where the 
intermediate Majorana neutrinos are far off mass-shell, the signal would be much weaker Other contributions such as the box 
diagrams etc in Ref. [5J are all of similar nature and much smaller than those considered here. 

IV. CONCLUSIONS 

We extended the existing literature to consider the AL - 2 rare decays of heavy mesons D, Ds, and B to a vector or pseu- 
doscalar meson final state. Since there have not been any direct experimental searches on these LV heavy meson decay modes, 
we calculated their decay branching fractions and proposed to search for them in the existing and forth coming experiments. 

We first re-evaluated the limits on the mixing parameters \Ve4\^, |V;j4p and |ye4V^4| from some decay modes which have 
experimental limits on the branching fractions for the heavy mesons, and obtained full agreement with those in Ref. |@]. We then 
translated the limits on mixing parameters and corresponding mass ranges to the relevant decays modes of D, D,, and B of our 
current interests as summarized in Table I. Finally, we calculated the decay widths and branching fractions for various LV decay 
modes by the limits on the mixing parameters and the heavy neutrino mass. We sampled the constraints on branching fractions 
as a function of the heavy neutrino mass as shown in the figures. 

Although the prevailing theoretical prejudice prefers Majorana neutrinos, the unambiguous signature to prove the Majorana 
nature of neutrinos is the experimental detection of a LV process. A detection in one of the LV heavy meson decay modes studied 
in our analysis would imply LV and hence the existence of a Majorana neutrino. At present, about 2.4 x 10^ D^D 1, 17,1 events 
and 5.5 X 10^ Dl^'D'^ jlSIl events have been collected by CLEO collaboration. So these decay modes D, — > eep, D, ejxp and 
D efip et al. which we studied in the Table I might show up in the current experiments if the mass of the heavy neutrino is in 
the range 1 GeV < m4 < 2 GeV. But those B decay modes for the range of the heavy neutrino is 2 GeV < m4 < 4 GeV cannot 
be detected presently due to small branching fraction. Fortunately, in the forthcoming LHC experiments, around 10" ~ 10'^ 
B meson events are expected which will provide us with chances of discovering all the LV B decays modes we studied. 
Therefore, we may have the opportunity to discover the LV process of heavy mesons B, D and via the distinctive channels of 
like-sign dilepton production with no missing energy. Hadron colliders may serve as the discovery machine for the mysterious 
Majorana neutrinos. 
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